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Abstract 

We estimate the Baryon Asymmetry of the Universe (BAU) produced in an inverse seesaw model containing extra 
light singlets, and with lepton number conservation prior to the electroweak phase transition. A CP asymmetry 
6 ~ ^(1) i s required to obtain a large enough BAU. We discuss the relation between the baryon and WIMP relic 
£N| . densities in baryogenesis scenarios using the out-of-equilibrium decay of a baryon-parent of mass M : when baryon 
£jq| number violation freezes out, the remaining density of baryon-parents is ~ M/m^x the WIMP relic density. So the 
^ baryon/WIMP ratio is ~ eM/mw- A natural explanation of the similar WIMP and baryon densities could be that 
■ CP violation is of order the ratio mw /M. 
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^ 1 Introduction 

In recent years, there has been some interest in relating the cosmological number density of baryons[l, 2, 3], with the 
\ number density of dark matter particles[4]. A popular approach, refered to as Asymmetric Dark Matter (ADM) [5, 6, 7], 
is to implement a single mechanism that generates an excess of particles over anti-particles, which materialises both in 
the baryons and in the dark matter. However, it is not necessary to relate the amounts of Dark Matter and baryons, 
1 but rather, to explain that both have their observed values. This (subtle) distinction is relevant, because the relic 
density of a weakly interacting particle (WIMP) is naturally of order the observed dark matter density. So it is 
■ not clear what is gained by dropping the "WIMP miracle" , to link the dark matter number density to the baryon 
\ asymmetry, which depends on an arbitrary CP violation parameter. Scenarios which implement the link in the inverse 
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direction [8, 9], explaining the baryon number density from the dark matter density, could be more interesting. But 
they require a CP asymmetry that is naturally 0(1). The aim of this paper is to follow the "subtle" distinction 
mentioned above: we assume that the Dark Matter are WIMPs (so their observed relic abundance is already natural 
and requires no explanation), and attempt to build a baryogenesis model, which naturally generates, in a similar 
■ cosmology, the observed baryon asymmetry. 

We consider an "inverse seesaw"-like model [10], extending the Standard Model(SM) with TeV-scale electroweak 
singlet Dirac neutrinos ip, additional light singlet fermions s and a L = 2 singlet scalar 4>. The ip participate in Yukawa 
interactions with SM neutrinos, and the </> have (L conserving) interactactions with ip and s. For generic choices of 
couplings, lepton number violation arises spontaneously when <fi gets a vev, after the electroweak phase transition, 
giving Majorana masses to the doublet SM neutrinos. However, there is a particular limit where lepton number is 
conserved and the SM neutrinos are Dirac, whose baryogenesis prospects are explored in section 3.5. Section 2 reviews 
our model and notation. Since we are building a model, we would like it to exhibit as many interesting features as 
possible. Our model can fit the observed neutrino mass differences. However, the additional light singlets cannot fit 
the reactor neutrino anomaly[ll] and Big Bang Nucleosynthesis; this is briefly discussed in section 2. 

The Baryon Asymmetry of the Universe (BAU) is generated via leptogenesis [12] at the TeV scale. As the 
temperature of the Universe drops below the mass of the singlet ips, their number density is depleted by annihilations 
and decays. We require this to occur prior to the electroweak phase transition. They can decay to SM neutrinos 
v and a Higgs, or to light singlets s and cf>, and CP asymmetries can arise in these decays. Since lepton number is 
conserved until <f> gets a vev, these lepton asymmetries would be of equal magnitude and opposite sign in the doublets 
and singlets. The asymmetries can survive in the plasma, once "washout interactions", exchanging lepton number 
between the doublets and singlets, are out of equilibrium. Section 3.3 estimates that the co-moving number density of 
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ips remaining, when this occurs, is "naturally" of order that of WIMPs. So to "naturally" obtain similar relic densities 
of baryons and WIMPs would require to "naturally" obtain 0(1) CP asymmetries in the tjj decays. As reviewed in 
section 3.2, large CP asymetries can be arranged by taking the tp masses to be quasi-degenerate[13]. The lepton 
asymmetry in Standard Model neutrinos then can be partially transformed to baryons by non-perturbative SM B+L 
violation. Dark matter is assumed to be some other WIMP, with a "usual" relic abundance. So this scenario gives 
similar abundances of WIMPs and baryon- parents, but similar dark matter and baryon number densities only arise if 
there is an 0(1) CP asymmetry in the parent decays. Unfortunately, such a large CP asymmetry is not easy to obtain 
in our model. We comment on this scenarios ability to relate the baryon and dark matter densities in section 4. 



2 Review, Notation and Masses 
2.1 Observations 

The mass density of baryons B in the Universe today, as extracted from a canonical AC DM model [14], is 

mJriB — n-s) 9 (iig — n-s) , 

pV Si = il B h 2 = 0.02255 ±0.0054 => y B = — — 8.7 x 10~ n (1) 

p c s 

where s = 2w 2 g*T 3 /45 is the entropy density, p c = 3Hq / (8irG) is the critical density to obtain a spatially flat Universe, 
and Ho — 74.2 ±3.6 km/sec/Mpc = lOO/i km/sec/Mpc. The cold dark matter relic abundance, for the same canonical 
ACDM cosmology with three light neutrinos, is 

^DAih 2 = 0.1126± 0.0036 (2) 

so that the ratio of dark matter particles to baryons is 

XEH^ 5 JIhL- , (3) 

where m p is the proton mass. Since m p mostly arises from QCD, and the dark matter should not have strong 
interactions, we assume that m p and to_dm are unrelated, and focus on obtaining similar baryon and dark matter 
number densities today: Yb/30 ^ Ydm & 30 Yg. This definition of "similar" allows a dark matter candidate in the 
mass range accessible to most direct detection experiments. 

Notice that the ratio (3) could be different in more compliciated cosmological models, such as those with extra 
sterile neutrinos. The CMB data alone allows the dark matter density to increase with the number of rclativistic 
species present at recombination (counted as the number N e ff of neutrino flavours): SQdm ~ A(SN e f fflDM /N e f /) 
[14]. Including other observations, such as H , favours N eff = 4.34^0^(68% CL) [14, 15]. An increase in Qdm 
would usefully increase the ratio (3), various laboratory neutrino anomalies [11] could benefit from light neutrinos 
other than the three of the Standard Model, and this paper constructs a neutrino mass model with additional light 
singlets. However, cosmological data also constrains the light neutrino mass scale. An additional cV sterile neutrino 
that could address the reactor neutrino anomaly should respect cosmological bounds on hot dark matter. It was 
shown by [15] that such a eV-sterile is cosmologically consistent if there are additional lighter neutrino species to delay 
matter radiation equality, and an 0(0.01) asymmetry in the v e density to allow Big Bang Nucleosynthesis (BBN) 
[16, 17] to produce the observed relic abundances of light elements. Our model allows several sterile species, but the 
baryon asymmetry is related to the asymmetry in SM neutrinos, so we do not obtain the large v e chemical potential 
required by [15, 17]. Furthermore, if our extra steriles decouple prior to the Electroweak Phase Transition (EPT), 
their contribution to the radiation density at BBN is suppressed by a factor g*(TsBN)/ ' 9*(Tept), where g*{T) is the 
number of Standard Model degrees of freedom in the plasma at T. 



2.2 Models 

An explanation for the ratio (3), which has benefited from recent interest [6, 7] could be that a single asymmetry 
controls both the baryon and dark matter relic number densities. If this ratio, given in eqn (3), should be exactly 
one, then a comparatively light cold dark matter candidate is required. A model where the DM gets a mass after the 
EPT is [18]. The observed difference in mass densitites could also be due to inequivalent cosmological histories after 
the formation of the asymmetry. Some early discussions of asymmetries in the dark matter density are [5, 19]; for a 
review of recent Asymmetric Dark Matter (ADM) scenarios, see e.g. [20]. Asymmetric leptogenesis models [7], and 
links to Dark Energy and inflation [21, 22] also have been explored. ADM models explain why the dark matter and 
baryon number densities are similar, but not why they have have the observed value. 
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An alternative approach [8, 9], attaches importance to the " WIMP miracle" (see eg [4]): the relic abundance of a 
massive particle, which annihilates with its anti-particle via a weak cross-section, gives SloMh 2 of order the observed 
value. To explain the ratio (3), one then only needs a natural scenario giving a relic number density of baryons minus 
anti-baryons, which is similar to the relic number density of WIMPS plus anti- WIMPS. It is well-known [1] that 
to produce an asymmetry, CP violation and non-equilibrium are required, whereas only a departure from thermal 
equilibrium is required for the particle-antiparticlc-symmetric relic density of WIMPs. So if the non-equilibrium 
is similar in the baryon and dark matter production, and the CP asymmetry for baryogenesis is one, then similar 
number densities could be expected. For instance, in [8], the baryon asymmetry is produced in the out-of equilibrium 
annihilations of the dark matter particle, which works for selected mass ranges of the particles involved. In this 
paper, the lepton asymmetry is produced in decays. The parent particle of the asymmetry first annihilates, then 
decays out-of-equilibrium, at the TeV-scale. This freezes-in the baryon-parent and dark matter densities at a similar 
temperature, so the departures from thermal equilibrium in both processes may be similar. Indeed, the co-moving 
density of a weak-scale unstable particle, when inverse decays freeze out, is similar to the relic density of a weak-scale 
particle that annihilates. So an 0(1) CP asymmetry is required to get similar baryon and dark matter densities today. 

We construct a TeV-scale model of out-of-equilibrium-decay leptogenesis [12]. Various mechanisms have been 
studied for leptogenesis 1 in TeV-scale seesaw models [25], such as degenerate decaying singlets [13], or adding extra 
particles [26]. Heavy singlet fermions in the "inverse seesaw" pattern have been studied by various people [27, 28, 
29, 30]. Our model differs from the usual inverse seesaw in that we have additional light singlets[31], so we can 
generate a non-zero asymmetry in SM lcptons, despite that the model conserves lepton number. Our leptogenesis 
scenario therefore differs from that of [27], who also study a lepton number conserving scenario, but without extra 
light singlets. They obtain a sufficient baryon asymmetry by turning off the sphalerons before inverse decays go out 
of equilibrium, when various asymmetries are present and changing. Here, we explore whether the baryon asymmetry 
can be obtained in a more "dynamics-independent" way, due to the presence of extra singlets. The Lagrangians of [28] 
contain "hard" lepton number breaking, so the relevant interactions for leptogenesis are somewhat different, because 
a net lepton asymmetry in generated. An advantage of the hard breaking is that it can naturally provide the small 
mass splitting between the heavy decaying singlets [28]. The model of [29] generates the asymmetry in the TeV-scale 
singlets. We focus on a lepton number conserving model, because it is for such models that a "hidden sector" (as 
provided by the extra light singlets) is useful for storing an asymmetry [32]. Futhermore, the natural scale of lepton 
number violation in the inverse seesaw is below the electroweak scale, so if lepton number violation is spontaneous, 
one could expect it to arise after the electroweak phase transition. 

The inverse seesaw [10] has received attention because it generates small neutrino masses due to new particles with 
TeV-scale masses and couplings of a "natural" size. Such new particles may therefore be kinematically accessible to 
the LHC [33], or induce detectable [34] non-unitarity and/or Lepton Flavour Violation. We wish to explore, in models 
where the small lepton number violation is spontaneous, the prospects of additional light singlets and leptogenesis. 
A careful study [27] (without additional light singlets s), showed that the baryon asymmetry can be generated by 
turning off the sphalerons during the decays. We had hoped that baryogenesis would be easier with additional light 
singlets, but the parameter space where we estimate that it could work is not large enough to be convincing. 

Dark Matter has been studied in radiative inverse seesaw models [35, 36], where it can arise more naturally than 
in the inverse seesaw [37]. We recall that our model does not aim to provide a WIMP; we suppose that it arises in 
some other sector of the theory. 

2.3 Notation and Masses 

We consider a Lagrangian with several additional singlets, which can spontaneously violate lepton number. The SM 
doublet leptons I have Yukawa couplings to (two or three) " right-handed neutrinos" Nr, there is some larger number 
of singlet left-handed leptons {S, s}, and an L = —2 singlet scalar <p. The left-handed singlets who participate in a 
mass matrix M with the Nr are refered to as {S}, and the remaining ss can acquire masses once <f> develops a vev. 
Allowing all renormalisable interactions, the Lagrangian, written in two-component notation with all the fermions 
left-handed, is: 

C = C SM + {A* [IH]N C + M* N C S + y* 4>Ss + — </>SS + —<f>ss + —^N C N C + h.c.} + V{4>) , (4) 

H is the Higgs, [£H] — vH® — eH+, and N c , S and s have lepton numbers L = { — 1, 1, 1}. The potential V(<j>) causes 
(f> to develop a vev after the electroweak phase transition. We assume that such a potential can be constructed, and 
consider it no further. We also ignore the resulting majoron [38]. Generation indices are implicit; for simplicity, we 

1 Low temperature out-of-cquilibrium-decay baryogenesis scenarios have also been considered; the analysis of [24] was useful to us. 
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suppose two generations of N, S and s, so the matrices A, y, Y, X, Z and M are two by two. We denote the eigenvalues 
of A as Ai, A 2 (and similarly for other matrices). We will see that Ai <~ A2, y\ ~ yi and M\ <~ M 2 to obtain a large 
enough CP asymmetry for leptogenesis. 

Prior to the Electroweak Phase Transition, there are two gauge singlet Dirac fermions, ipi 7 ip2'- 



ipi = 



Si 



(5) 



whose masses Mr ~ TeV are the eigenvalues of M. The Mj will be taken degenerate to obtain a large enough baryon 
asymmetry, so when the mass difference is irrelevant, we write M$. 

In the presence of vacuum expectation values v = (H) = 174 GeV and u — ((f) ^ v, the neutral Majorana mass 
matrix can be written, in two component notation, as 



~{v L N c Ss) 
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h.c. 



(6) 



where mu = Xv, fix = Xu, fi v = yu, fiz = Zu, and fiy = Yu are two by two matrices. 

We are interested in the limit where the eigenvalues of M are much larger than the other entries in the mass matrix. 
In the case of the usual inverse seesaw, which contains no light singlets s, the determinant implies that the light active 
neutrino mass matrix is <~ mrjM -1 fiYM T ~ x m?Q. In the case where extra singlets are present, the determinant of 
the mass matrix in (6), in one generation, is m 2 D (fiYfix — My)- We neglect from now on the coupling Z, because its 
contributions to the mass matrix are unimportant, and its effects in the baryogenesis scenario are similar to those of 
Y, which will be unimportant. We will see in section 3.3 that the couplings in y and/or A must be small 2 (£3 10~ 5 ), to 
prevent the washout of the lepton asymmetry. We focus on two simple limiting cases: yi -C Yi,Xi ^ 1 , which gives 
Majorana masses for the active neutrinos, and the case Yi = Xi = 0, which gives Dirac masses between the vl and s. 

In the Majorana case with yi <IC Yi,Xi, the active neutrinos have a usual- inverse-seesaw- like mass matrix m v ~ 
rri£iM"VyM T " 1 mJ, whose eigenvalues satisfy 



^atm^sol 

eV 2 



= 0.3 



A1A2 



(3TW)_ 2 V Y.Y^f 
M X M 2 



(7) 



<j> should have a mass and vev well below the Higgs vev v = 11 A GeV, to ensure that lepton number is conserved prior 
to the electroweak phase transition. However, we will be forced to take the vev of order v (to allow A, small enough 
to give a sufficient abundance of ips after freezeout of washout interactions). We will also require Ai ~ A 2 , to obtain 
a large enough CP asymmetry, so we will suppose a mild hierarchy in Y to generate the atmospheric - solar splitting. 

The light steriles s have a mass matrix ~ \xx — X {</>). In the one generation case, s has a mixing angle with 
the active neutrino of order ^ y m£i/(M^x), which allows the decay s — > iv at tree level. We focus on parameters 3 
X ^ Y, such that the singlets are innocuous: their masses are ^ GeV, they annihilate efficiently (to majorons), and 
otherwise can decay. Since the ss decouple prior to the Electroweak Phase Transition, their temperature at T ~ fix 
is suppressed with respect to the photons, so they should not over-contribute to the radiation density. 

The active neutrinos share Dirac masses ~ mofiy/M with the singlets s, \iX = Y = Z = Q\n. the Lagrangian (4). 
We consider this limit only in section 3.5, because it allows more parameter space for our leptogenesis scenario 4 . By 
assigning L = to <fi, and -1 to s, lepton number clearly is conserved also after both (f> and H get vevs. However, since 
<j> was introduced to spontaneously break lepton number, it is peculiar to not give it lepton number. We nonetheless 
assume that (j) does not mix with H, and do not discuss potential constraints from Higgs physics. 



3 Thermal History 

One of the roles of ip is to generate the Baryon Asymmetry of the Universe (BAU). In its decays, it can generate an 
active versus sterile lepton asymmetry and the excess of active leptons will be partially transformed to baryons by 

2 We therefore do not expect bounds on our model from lepton flavour violating processes such as fi — ¥ e-y, or non-unitarity [34], 
3 With X ~ 10~ 10 , we can obtain a singlet s with eV mass and ~ 0.1 — 0.01 mixing angle, that could fit the reactor neutrino 

anomaly. However, we do not explore this parameter space, because out model cannot give the large lepton asymmetry required to make 

it cosmologically acceptable. 

4 Mixed scenarios can be envisaged, for instance where the atmospheric mass is Dirac, and the solar mass is Majorana. We do not study 
this tuned example because we estimate that out-of-equilibrium-decay-leptogenesis does not occur. 
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sphalerons [23]. This is discussed in section 3.2. The asymmetries produced in tp decay can survive in the plasma 
once "washout interactions" (such as inverse decays £H — > tp and (ps — > tp, and scattering IH — > (ps) go out of 
equilibrium. This "freeze-out" temperature Tbau is estimated in section 3.3. The baryon asymmetry present today 
will be generated in the decays of the tps remaining at Tbau- 

Decays and inverse decays are not the only interactions which change the tp number density; tp + tp can also 
annihilate, via the coupling Y to <f> + (p. The annihilation rate is faster than the decay rate at temperatures just below 
the tp mass, so we start by discussing annihilations below in section 3.1. Annihilations and inverse decays will freeze 
out at similar temperatures. 

3.1 iftip — > (f)(f) annihilations 

At temperatures T <; Mj, an equilibrium abundance of tpi and tpj will be present, because they can be produced via 
Yukawa interactions involving A (for instance in tRth — > H* — > £tp, at a rate T ~ h% A 2 T/(47r), which is fast compared 
to the Universe expansion H). 

As the temperature T drops below their mass, the singlets can annihilate via their Y couplings, or decay to a light 
lcpton and a scalar via their A or y couplings: 

tp I +^j-^ ( p + ^ (8) 
il)i-*l + H ,8 + ~$ (9) 

We focus first on annihilations, because we envisage 1 <~ » yj, Afe ^ 10~ 4 , which implies that initially, at T £ Mj, 
the tpis are more likely to annihilate than decay. 

The non-relativistic annihilation cross-section in the centre of mass is 

™^ + ^^) = Sf ^ 

where all the singlet masses are approximated as M^. We neglect the Sommerfeld enhancement due to (p exchange 5 
because the tps are almost relativistic; the enhancement can be estimated [39] to be a few percent. The number 
density of tpis remaining when annihilations freeze out, and the freeze-out temperature T ann = M^/z ~ M^,/20 can 
be determined from 

m IT \^n,^TJf r T \ ~. hi 'V ' 9*{Tann)Tl nn . . 

n^{±annpV ^ H(I ann ) ~ (llj 

where g* is the number of degrees of freedom in the plasma, and g*(TeV) ~ 100 in our model. 

3.2 Decays and the CP asymmetry 

After the tp annihilations freeze out, the number density of ^s and tp will continue to drop, because they decay. The 
model conserves lepton number, so the total lepton asymmetry produced in these decays is zero. However, since the tps 
decay both to SM particles H, £, or to singlets s, (p, equal but opposite lepton asymmetries in the singlet and doublet 
sectors could be produced. To obtain a baryon to dark matter ratio that is 0(1), these CP asymmetries will need to 
be large, so we explore this limit below. 

The tpjs, which carry L = 1, can decay at tree level to an l a and an H, or to a (p (L — 2) and a Hp, at partial rates 

t a + H) = , rty I ^ap + t) = ^M I (12) 

See the first diagrams of figures 1 and 2. These decays are fast, compared to the expansion rate H, for T <, 
{|A|, \y\} ,y M i nipi/ 20, so at all temperatures relevant to us. 
CP violating asymmetries for a given final state: 

(i,/) = n^i -> tgH) - rgj l a H) (J , a) = rjtpr -> s a -$) - rg 7 -> Sa <b) 
r(^/ -> a/0 + r(^j -> ail) ' e " " r(v>/ -> a//) + r(^j -> o/i) 

can be obtained from the interference of tree with loop diagrams, if the couplings are complex and some particles 
in the loop can be on-shcll. In figures 1 and 2, are drawn the subset of one- loop diagrams which we include in our 



5 We thank Marco Circlli for discussions of this issue. 
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Figure 1: diagrams contributing to the CP asymmetry with final state £ a + H . 




Figure 2: diagrams contributing to the CP asymmetry with final state s a + <j>. 



calculation. We are interested in e — > 1, which can be obtained in the decays of quasi-degenerate singlets from the 
illustrated wave-function renormalisation diagrams[13], so we neglect the additional vertex corrections diagrams. 

We neglect the effects of X and Y couplings in the decay rate and CP asymmetries, despite the hierarchy in 
magnitudes: X, Y <~ 1, y, A ~ 10 -4 . This is because decays involving X and Y are suppressed by three-body final 
state phase space, and we were unable to find significant effects of X and Y in the CP asymmetries. 

In the decay of tpi, the total asymmetry in the doublets, summed on lepton flavours a, is 

Aid) _ y Aid) _ y Im {^iKjVpiV}A (u] 
^ » ^ SrrpAtA + yty],, l- Xj { > 

where xj — Mj/Mf. This contribution arises from the interference of the tree decay with the last loop of figure 1. 
Similarly, the total asymmetry in the singlets is 

M,s) _ V A 1 -*) = V" Im {v**iy>*ApA*pi\ V^l , , 



In both the doublet and singlet sectors, there can be asymmetries in the individual flavours a, which arise from 
the middle diagram of figures 1 and 2. We neglect these contributions, and focus on the flavour-summed singlet and 
doublet asymmetries of eqns (14) and (15), because its the total (=flavour-summed) doublet lepton asymmetry which 
the sphalerons transform to baryons 6 . 

The asymmetries given are for the decays of any of the ipi. To obtain a large enough CP asymmetry, at least two 
of the il>i must be very degenerate :M/ — Mj ~ T/, so the CP asymmetries in the decays of at least two of the ipis will 
contribute to the baryon asymmetry. It is therefore fortunate that the CP asymmetries do not vanish when summed 
on /: in a two generation model, in the limit where Mi —> M 2 , the CP asymmetries are equal e^ 1 '^ = e^ 2 '^, rather 
than opposite. 

As anticipated, the sum of e^ 1 '^ and e^ J ' s ' vanishes because the model conserves lepton number. It will therefore 
be important that no interaction be fast enough to equilibrate the singlet and doublet sectors, as the asymmetries are 
produced and until the Electroweak Phase Transition. This is discussed in section 3.3. 

It is helpful to be able to estimate the size of these CP asymmetries. This scenario of out-of-equilibrium decay 
leptogenesis at the electroweak scale requires a large CP asymmetry e <; 0.1. So from eqn (14), all the elements [A]/ Q 
and [y]ip must be comparable, for all a, I, (3, and the degenerate :M\ — M 2 <; T. For phases to maximise the 
asymmetry, large mixing angles, and comparable eigenvalues Ai <~ A 2 , yi ~ yi, 



e 



y(X{ M I ^ 1 



< 77^ (16) 



16tt y\ + 2A 2 Mj - Mj ~ 4D 
where the last approximation expresses the mass splitting in units of Tj 

Mi - Mj = DYi (17) 



6 The flavour asymmetries can nonetheless generate the observed baryon asymmetry, due to time-dependent effects [27] 
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and the coupling constant combination was maximised by taking 

2/i~\/2Ai . (18) 

Recall that the formulae (14) and (15) are valid for Mj — Mj > Tj, so e ~ 0.1 is barely consistent. Since the mass M is 
the only interaction linking the singlet and doublet sector, loop corrections should not destabilise the small splitting. 

3.3 Washout of the asymmetry by inverse decays and scattering 

As the decay, they produce an asymmetry in doublet leptons. This asymmetry will survive, and be redistributed 
throughout the plasma by the interactions in chemical equilbrium [40], if all the interactions which can destroy this 
asymmetry arc out of equilbrium. Such destructive interactions are refered to as "washout interactions", and in 
our model, are the inverse decays £H — > tp and <ps — > ip, which transfer lepton number back to the tps, and the 
Hi — > (j)s scattering (see figure 3) which exchanges lepton number between the doublets and light singlets. We follow 
the tradition 7 of treating the on-shell-^ part part of figure 3 as inverse decays, and including the remaining " Real 
Intermediate State -subtracted" part as scattering far below the tp pole. 




Figure 3: Dangerous scattering process I + H — > s + <fi 

For making estimates, it is convenient to have thermally averaged rates, as opposed to the reaction densities 8 used 
in Boltzmann Equations. So we define the thermally averaged inverse decay rates as 

n e (l ) ne{l ) 

where ng(T) ~ ggT 3 /ir 2 7 the tp decay rates are given in eqn (12), and the equilibrium abundance of a non-relativistic 
particle X is 



m x T 
2tt 



3/2 



e- mx ' T (20) 



where gx is the number of spin degrees of freedom of X (notice that our number densities describe particles, not 
the particles + anti-particles). Comparing the inverse decay rate of eqn (19) and the annihilation rate of (11) to the 
Hubble expansion, shows that they freeze out at a similar temperature 

T BAU * | (21) 

The annihilations, discussed in section 3.1, will freeze out first for 

M 3 ! 

\[YY^ n \ 2 <^Ay ? A 2 AAKA 2 ■ (22) 

We will see that this condition is satisfied, so annihilations stop depleting the ip number density before the inverse 
decays allow a lepton asymmetry to survive. 

The number density of ips remaining when both inverse decays (T{£H — > tpj)) and (T(<fis — > ipi)) arc out of 
equilibrium is 

(rp N • / H(TBAU) (rr \ H(T B AU) f rp x\ / oq \ 



7 See [3] for a discussion, and references, of this curious separation. 
8 see, for instance, [3] for definitions. 
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where T and H are given in equations (12) and (11). 

To protect a doublet lepton asymmetry, the inverse decays from the singlet and doublet sectors must both be out 
of equilibrium. To see this, recall that the heavy Dirac ips have L = 1, and there are equal numbers of ips and ips to 
begin. As they decay, opposite lepton asymmetries develop in the doublet leptons and light singlet sector. Consider 
now the case that inverse decays from the doublets are out of equilibrium, but that the asymmetry in light singlets 
can be transfered back to the heavy Dirac leptons by (ps — > ipi . Then the anti-asymmetry from the light singlet sector 
is transfered to the ips and ips, who transmit it to the doublets via decays, which are always in equilibrium. 

Away from the ip mass pole, the scattering cross-section for I + H — > s + <p is 

( . |y/3/A aJ | 2 Mf \ypiX aI \ 2 

where, after the arrow, the internal line momentum is neglected. This is an acceptable approximation at temperatures 
T & M//20, where the inverse decays are out of equilibrium. With this approximation, the thermally averaged 
scattering rate T ~ erne is out of equilibrium for 

| m A Q/ | 2 ^4xlO- 11 ^. . (25) 
3.4 Putting it all together in the Majorana case 

The "washout" interactions go out of equilibrium at Tbau (given in eqn (21)), after which an asymmetry in the SM 
leptons is produced. This should occur prior to the electroweak phase transition, which we take to occur at Tept ~ 
150 GeV [41]. Requiring T B au > T EPTl imposes a lower bound on the singlet mass scale 

M £ 3 TeV . (26) 

There is also an upper bound on M, as a function of A and Y, from requiring that the inverse seesaw give the observed 
neutrino mass differences. This condition was given in eqn (7). Since a large enough CP asymmetry requires large 
mixing angles and comparable eigenvalues in A and y, we can restrict to the one generation case: 

„,„„= „v^Wp^>yw> eV m 



lcrV V M 2 J 

And finally, the requirement that ipip — > (pep annihilations freeze-out before the inverse decays £ a H — > ip, imposes the 
relation between Y and A given in eqn (22). Expressed in terms of eigenvalues, we obtain 

# * 10_S (lF*) • (28) 

Electroweak sphalerons will partially transform the B — L asymmetry produced in SM fermions, Yb-l, to a baryon 
asymmetry[40] Yb = (12/37)Yb_£. So the baryon asymmetry produced in our model will be of order 

^ U nt-nj „ 12 4MTW) „ 10 _M±_ 10^1 . , 

B "37 s "37 s e - 4xiU 3TcV|A 2 | 2 i? ' ^ 

where 4n^, = Ej(n^ )J +n^ I ), and we used eqns (23) and (16). It is marginally possible to obtain the atmospheric mass 
difference and the observed baryon asymmetry, with M = 3 TeV, Y<i saturating eqn (28), and (</>) ^ v. In figure 4, the 
available parameter space is plotted as a function |j/x | 2 / 10 9 in eqn (17). 

3.5 Dirac limit 

As seen in the section above, for Majorana light neutrino masses, the parameter space of our model is very restricted. 
In this section, we explore the Dirac limit of the light neutrino mass matrix (6), where all Majorana mass terms are 
zero, /ia — ^ ,ie {A, Y, Z}. Without Majorana mass terms one can redefine the lepton number of {N c ,S,s,(p} 
to be { — 1, 1, —1,0} so that also after the phase transition of <f> lepton number is conserved. Neutrinos become Dirac 
particles with si being their right handed Dirac partners. The light mass eigenvalue is 

< = ^ (30) 
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Figure 4: Available parameter space generating the observed baryon asymmetry and neutrino masses, of Majorana 
type, in the inverse seesaw model with extra light singlets, described by the Lagrangian (4). 

which is only suppressed by jj compared to -A^ in the Majorana case. This will significantly enhance our allowed 
parameter space. But it is also clear that hirarchy in neutrino masses must now come from a hierarchy in the Yukawa 
eigenvalues Ai < A2, and/or the same for y. 

Before looking at the parameter space in the Dirac model, we shall look at its thermal history, comparing to the 
Majorana case. The first thing to realize is, that annihilation processes do not exist. The only processes that washout 
CP asymmetry are scattering, as in fig. (3) and inverse decays which, like in the Majorana case, determine the moment 
of efficient CP violation. We assume the inverse decay rates from the various flavours are comparable, despite the 
mild hierarchy in A and/or y. After the inverse decays become slow, all heavy parent particles will decay and produce 
a CP asymmetry that is not washed-out, so that the first part of equation (29) remains valid. Therefore the thermal 
history in the Dirac limit does not differ significantly from the Majorana case. The value of e, however, will be reduced, 
because of the Yukawa hierarchy. To be concrete, in figure 5, we take e = (compare to the Majorana case, eqn 
(16)). 

One should also have a careful look at Big Bang Nucleosynthesis because it restricts the number of thermalised 
neutrinos at T ~ MeV: N v ^ 4. It was shown in [42] that right handed Dirac neutrinos are not in equilibrium at 
T ~ MeV because their interaction rates are suppressed by a factor (^f) compared to that of left handed neutrinos. 
Consequently the number of neutrinos in thermal equilibrium at BBN is not changed by our model with respect to 
the Standard Model. 

Let us now turn to the parameter space of the Dirac limit. In the Majorana case the allowed parameter space is 
very small because asymmetry production favours small y, in order to produce a big density of ip when inverse decays 
freeze out which can then be transformed into enough baryon asymmetry, even if e < 1. But neutrino masses prefer 
big y couplings, because the neutrino mass term is already suppressed by 1/M 2 . This forces us to take quite large 
values for u, the vev of <f> in order to enhance neutrino masses without limiting asymmetry production. In the Dirac 
limit neutrino masses are less suppressed, so that the conflict between baryon asymmetry and neutrino masses is less 
critical. We can realize our model with moderate values for u « 6GeV, far away from electroweak phase transition, 
while parent masses are at M = 3TeV and do not have to be more degenerate than Mj — Mj ss 5T(tpj — > all). These 
values were obtained for a ideal constellation of our parameters, i.e. all phases in our Yukawa couplings A and y are 
of order one and that 2A| ~ y\. However in this Dirac limit small variations from these ideal conditions are possible 
as the parameter space is big enough to compensate for non ideal conditions. A plot of possible parameter space is 
given in fig. (5). 

4 Discussion 

This paper considered baryogenesis at the electroweak scale, in an out-of-equilibrium decay scenario. This was im- 
plemented in an inverse seesaw model with extra light singlets s, described by the Lagrangian (4). The aim was to 
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Figure 5: Constraints from the Baryon asymmetry and neutrino masses on y\ as a function of the parent particle 
mass scale M$ give the allowed parameter space in the Dirac limit. The upper bound from Baryon asymmetry (solid 
line) varies with the degeneracy Mj — Mj w DT(tp] —> all). For less degenerate parent particles, the slope reduces 
and the allowed parameter space shrinks. The dashed line is the lower bound from neutrino masses and varies with 
the vacuum expectation value u. For smaller u, the slope rises and the parameter space shrinks. 



relate the baryon asymmetry to the dark matter abundance, and possibly fit various neutrino anomalies with the extra 
singlets. Similar Dark Matter and baryon abundances were supposed to arise, because, if the Dark Matter is made of 
WIMPs, and the baryon asymmetry produced in electroweak-scale out-of-equilibrium decay, then both relic densities 
are controlled by electroweak-scale interactions going out of equilibrium. The baryon asymmetry is also proportional 
to a CP asymmetry e, see eqn (29). 

The scenario for generating the baryon asymmetry in our model is outlined at the end of the introduction. It does 
not work very well. The first difficulty is that it does not naturally give the large CP asymmetry that is required, 
contrary to models with explicit lepton number violation (reviewed in section 2.2). More importantly, the allowed 
parameter space where it produces a sufficient baryon asymmetry and correct neutrino masses is tiny (see figure 4), 
and unattractive. It requires a small mass splitting between the heavy singlets M% — M\ ~ T, but the formula for 
e, eqn(14), becomes unreliable in this limit. Also, the inverse seesaw is attractive because it allows new LHC-scale 
singlet fermions with large Yukawa couplings, and spontaneous lepton number violation due to a small vev. However, 
our model requires small yukawas A and y, to ensure that interactions which wash out the baryon asymmetry freeze 
out soon enough. This forces a large lepton number violating vev, of order the Higgs vev. The masses of the extra 
light singlets s are arbitrary, but would be ^ GeV, unless a steep hierarchy is tuned into the couplings (see discussion 
after eqn (7)). A larger successful parameter space is obtained (see figure 5), if the Majorana couplings of the inverse 
seesaw model are set to zero. Then the active neutrinos obtain Dirac masses with the light steriles. However, then 
the motivation for the singlet scalar (whose vev spontaneouly broke lepton number in the Majorana case) , is no longer 
clear. Since it no longer carries lepton number, it could mix with the Standard Model Higgs, and some analysis would 
be required to determine its experimental signatures. 

Despite the flaws in our model, it could be unsurprising, if not natural, to obtain similar baryon and dark matter 
relic densities, when the dark matter are WIMPs, and the baryon asymmetry is produced in out-of-equilibrium decay. 
In such scenarios, the baryon asymmetry is produced after the freeze-out of washout interactions. This condition can 
be roughly estimated as 

T(all -> ip) ~ T(ip ->• all)e~ M/T < H(T BAU ) (31) 

where tp is the decaying parent of the BAU, of mass M, and H is the Hubble expansion rate. Notice that inverse 
decays are Boltzmann-suppressed, because the light decay products have difficulty to find, in the thermal bath, the 
energy to produce a heavy ip. 

This Boltzmann suppression is reminiscent of the freezeout of annihilations of a WIMP which can be estimated 
to occur when 

n x (T DM )(av) £ H(T DM ) (32) 
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The e m l T appears here in equilibrium number density of a non-relativistic particle, see eqn (20). 

From these two estimates, it is straightforward to estimate the ratio of co-moving densities 9 of BAU parents tjj to 
WIMPs x as 

Y± 1 / 9(Tdm) . . 

Y x ~ (2ttz) 3 / 2 A 2 m x y g(T BAU ) {66> 

where the decay rates and annihiation cross-section were normalised to the relevant particle masses T = X 2 M^,/(8tt), 
and (av) = g 4 /(8Trm x ), and we assume that the baryon-parents and WIMPs freeze out at comparable values of 
z = m/T ~ 25. At a given mass scale, this shows that the co-moving number density of BAU-parents can be easily of 
the same order as that of particles which annihilate. So weak-scale BAU parents naturally have the same abundance 
as weak-scale dark matter. 

However, all BAU parents do not neccessarily produce baryons. On average, a parent produces e baryons, where e 
parametrises CP violation, and is suppressed by a loop factor, mixing angles and possibly small couplings. This could 
suggest that e is naturally ^ 10 -3 . Eqn (33) predicts 

n B Yxj> _ n M^p 

<~ e— — Oe 

TlDM y x m x 

where the coefficient C could be 0(1), for A 2 <C g 4 - So to obtain )is ~ udm requires e <~ 1 for m x <~ or m x < 
for e < 1. It could be reasonable to suppose that the WIMP is an electroweak-scale particle, with gauge couplings 
g, and that the baryon-parent is a flavour-scale particle, with flavoured couplings A, in which case, such a ratio of 
couplings could be credible. Then for "natural" values of the CP asymmetry e 1CP 3 , and the flavour-scale M$ *J 100 
TcV, the baryon and WIMP number densities today are "naturally" comparable. 
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